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Summary. — We examine the prospects for the detection of Higgs mediated lep- 
ton flavor violation at LHC and at a photon collider in the minimal supersymmetric 
standard model with large lepton flavor violating mass insertions in the fi — r sec- 
tor constraining the parameter space with several experimental bounds. We find 
rates probably too small to be observed at future experiments if models have to 
accommodate for a neutralino relic density as measured by WMAP and explain the 
(g — 2) M anomaly: better prospects are found if these two constraints are applied 
only as upper bounds. The spin-independent neutralino-nucleon cross section in the 
studied constrained parameter space is just below the present CDMS limit while 
gamma rates from neutralino annihilation in the halo are strongly suppressed. 

PACS 11.30.Pb - Supersymmetry. 

PACS 14 . 80 . Da - Supersymmetric Higgs bosons. 

PACS 95.35.+d - Dark matter (stellar, interstellar, galactic, and cosmological) . 



1. Introduction 

The Higgs sector of the minimal supersymmetric standard model (MSSM) [1], es- 
pecially the heavy neutral Higgses A and H, play a prominent role in the physics of 
neutralino dark matter [2]. In some region of the supersymmetric (SUSY) parameter 
space neutralinos yield the desired amount of relic density by annihilating into fermions 
through the s-channel resonant exchange of neutral Higgs bosons h, H, A, the so called 
funnel region where vtia — 2m x ; besides, as dark matter is expected to be distributed 
in a halo surrounding our galaxy, neutralinos can scatter off nuclei in terrestrial detec- 
tors: the coherent scattering is mediated by scalar interactions through the s-channel 
exchange of squarks and i-channel exchange of the CP-even neutral Higgs bosons h and 
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H . These effects become sizable when squarks are heavy and tan (3 is large in reason of 
the enhanced Higgs bosons coupling to down-type fermions, especially for the b quark 
which has the largest Yukawa coupling receiving large radiative SUSY-QCD corrections 
at large tan [3. 

Once a source of lepton flavor violation (LFV) is present in the slepton mass matrix, 
for example the MSSM with the see saw mechanism for generation of small neutrino 
masses [3], non-holomorphic LFV Yukawa couplings of the type L R L J L H* are induced at 
loop level and become particularly sizable at large tan/3 giving rise to enhanced Higgs- 
mcdiated LFV effects [4] . The LFV mass insertions 5 l [ L — (rn 2 L y 3 /m 2 L , S RR = {m R ) K '/ m i?> 
where (m| R ) 13 are the off-diagonal flavor changing entries of the slepton mass matrix, 
are free parameters which allow for a model independent study of LFV signals. We 
introduce LFV in the model through the mass insertions dj? L RR — 0.5. This value 
ensures the largest rates in LFV processes and allow us to study the more optimistic 
scenarios for LFV detection; higher values contradict the mass insertion approximation 
as an expansion of propagators in this small parameters. 

Higgs mediated effects become interesting at large fi and tan/3 and low m^; further, 
if SUSY-QCD particles are heavy, Higgs effects are dominant also for neutralino dark 
matter. We thus consider the following real MSSM parameter space: 100 GcV < wia <1 
TeV, 20< tan/3 <60, 500 GeV< // <5 TeV, (the sign of \i is taken positive, as preferred 
by the SUSY explanation of the (g - 2) p anomaly), 150 GeV< M 1 ,M 2 <1.5 TeV (we 
do not impose any relation but let them vary independently allowing for gaugino non- 
universality at the weak scale), 1 TeV < M3 < 5 TeV (to have large masses for gluinos), 
1 TeV< mu 3 ,mD 3 ,mQ 3 <5 TeV (for the first and the second generation the soft masses 
are set to be equal, mjj i = m£> i = m,Q i = rriq, where i = 1,2 and m q is another free 
parameter which varies in the same range), 300 GeV< m^, rriE 3 < 2.5 TeV (for the first 
and the second generation the slepton soft masses are set to be equal, to^ = = rug, 

where i = 1, 2 and to;), -2< — — , — — < 2 while for first and second generation 

' £/' — mu 3 mr> 3 niE 3 — 

the trilinear scalar couplings are set to zero. 

We impose on the parameter space several experimental limits: (1) LEP, TEVATRON 
bounds on sparticle masses and the LEP bound on light Higgs; (2) present bounds 
on B-physics observables B(B — > X s j), B(B S — > ^i+/i~), Ato_b s , B(B — > tv); (3) the 
present experimental upper bounds on LFV processes B(j — > ^7), B(r — > fit]), B(r — > 
/i/i/i), (4) the exclusion limits on the neutralino- nucleus spin independent cross section 
from the CDMS experiment [6]; (5) the CDF exclusion exclusion limits in the (tua, 
tan fi) plane [7] . All of them are applied at the same time with the exception of the 
WMAP 3cr interval on relic density and {g — 2)^ anomaly for which we also relax the 
lower bounds: thus in the following figures the light gray (turquoise) points have only 
Vt x h 2 < 0.13 and aff SSM < 4 x 10~ 9 , the plus-shaped points dark-gray (indigo) points 
satisfies 0.09 < tt x h 2 < 0.13 and a^ SSM < 4 x 10~ 9 , finally, the squared points satisfy 
0.09 < VL x h 2 < 0.13 and 1 x 10~ 9 < aff SSM < 4 x 10~ 9 . For numerical computations we 
use the code DarkSusy [8] and the code FeynHiggs [9], while for the explicit formulas 
and further details on the experimental constraints we refer the reader to ref. [10]. 



2. Higgs mediated LFV at LHC and at a photon collider 



At high tan /3 the dominant production mechanisms for A, H at LHC is bb fusion due to 

m with Fe^ 

r(*^fc6) SM 



the TObtan/3 enhanced 66$ couplings. We calculate the cross section with FeynHiggs 
which uses the approximation a MSSM (bb -> *) = <j SM (bb -> $) r |!^"7 , where 
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Fig. 1. - Left: Scatter plot of the inclusive production cross section pp — > $ + X times the 
branching ratio of $ — > Tfi at LHC versus tua = A, H). Right: Scatter plot for the cross 
section of the process 77 — > Tfibb in photon-photon collision at y Js 11 = 600 GeV. The light 
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a SM (bb — > $) is the total SM cross section for production of Higgs boson with mass 
771$ via bb fusion: to obtain the value in the MSSM it is rescaled with the ratio of 
the decay width of the inverse process in the MSSM over the SM decay width [9], Wc 
calculate for each random model the product the cr{pp — s- $ + X) x £>(<(> — > Tfi). As 
masses and couplings of A and H are practically identical as discussed above, we have 
a(pp -t A + X) + a(pp -> H + X) ~ 2a(pp -> A + X). The scatter plot a(pp -> 
$ + X) x S($ — > t/.i) is shown in Fig. 1, left panel. We see that with the nominal 
integrated luminosity of 100 fb^ 1 per year models which satisfy both the relic density 
abundance and Aet M can give up to 10 events per year (squared points), up to 40 if we 
relax the condition on the lower limit of Aa^ (plus-shaped points) and up to 200-300 
relaxing both the lower limits (turquoise (light-gray) points). 

In 77 collisions the main production mechanism for $ = A, H is tt fusion while the 
bb is suppressed by a factor 3(l/3) 4 (m;,/777 r ) 2 ~ 0.1 which cannot be compensated by 
corrections to the b Yukawa coupling. In ref. [13] we studied in detail the fir fusion 
process 77 — > firbb where the Higgs boson is produced in the s-channel via a virtual 
fj,r pair and can be detected from its decay mode A — > bb. There we have shown that 
a good analytical approximation for the cross section is obtained using the equivalent 
particle approximation wherein the colliding real photons split respectively into r and fj, 
pairs with the subsequent /it fusion into the Higgs boson and that the effect of photons 
spectra can be neglected, we thus consider monochromatic photons with v /s 77 = 600 
GeV, and photon-photon luminosity 500 fb _1 yr^ 1 [13]. The scatter plot of the signal 
cross section versus 777,4 is shown in Fig. 1, right panel. Here the models which satisfy 
both the relic density abundance and Aa M (squared points) have maximal cross section 
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Fig. 2. - Left: Scatter plot for the spin-independent neutralino-nucleon cross section versus 
the neutralino mass. The area above the solid line is excluded the CDMS final results; the 
area above the dotted line is excluded by the 2008 CDMS search. The dashed and dot-dashed 
lines give the sensitivity reach of two phases of the XENON experiment. Right: Scatter plot 
for the particle physics factor entering the formula of the flux of gamma rays from neutralino 
annihilation in the halo, see the text. The light gray (turquoise) points have only Q x h 2 < 0.13 
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10 3 fb, which is too small. Relaxing the lower limits cross section values up to 2 x 10 2 
fb are possible, giving 10 events/year. 

3. — Neutralino dark matter direct and indirect detection 



The spin-independent neutralino-nucleon cross section in the limit of heavy squarks 

tan 2 j3 x Kf, where 



471-m^m 4 , 



and large tan (3 can be approximated as [11] a 

Nn and N13 are the lightest neutralino unitary mixing matrix elements, mjv the nucleon 
mass (neglecting the mass difference between the neutron and the proton) and Kf a 
factor which depends on nucleon form factors. The left panel of Fig. 2 shows the scatter 
plot for the spin-independent neutralino-nucleon cross section as a function of m x and the 
region excluded by CDMS [6] . We emphasize that CDF and CDMS limits are very mild 
constraints: the region excluded by CDF is practically excluded by the other constraints 
while the CDMS limit exclude only one plus-shaped point leaving untouched the regions 
preferred by WMAP and the (g - 2) M anomaly. The XENON100 experiment [12] should 
reach the sensitivity corresponding to the dashed gray (red) line in the Figure 1, left 
panel. Such sensitivity is able to cover the region with the highest cross section, m x > 300 
GeV, where there is large higgsino component. On the other hand the region preferred 
by (g — 2) M anomaly cannot be covered. We also report the prospected sensitivity goal of 
the XENON experiment with 1 ton detector mass [12], dot-dashed gray (red) line, which 
is 10 -11 — 10~ 10 pb for neutralino mass in the range 100 — 1000 GeV: practically all of 
the parameter space can be probed. 
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We also present the effect of our constrained parameter space on the flux of photons 
coming from neutralino annihilation in the halo which is a very active field in indirect 
dark matter detection. The flux of gammas expected from neutralino annihilation is 
generally given by F = $pp x $ astro where the second factor contains the astrophysical 
informations and is given by the integral of the squared of the dark matter density 
along the direction of observation, while $pp is the particle physics factor given by 
$pp(£ 7 > E th ) = (av}/(2m 2 x ) f™* dE(dN 1 /dE) where (av) is the thermal averaged 
cross section annihilation times the velocity of neutralinos and the dN^/dE is photon 
spectrum which is integrated over energies greater than E t h- Here we are interested 
mainly in the evaluation of the particle physics factor and in comparing it with other 
other studies in literature in the framework of the MSSM. We present $pp as a function 
of the neutralino mass in fig. 2, right panel, with threshold energy E t h — 1 GeV. For 
models with a relic density inside the WMAP interval the maximum value of the particle 
physics factor is 2 x 10~ 32 photons GeV" 2 cm 3 s" 1 for neutralino mass around 200 GeV: 
this value is two order of magnitude smaller than the typical values found in similar 
studies without the constraints from lepton flavor violation and the updated B physics 
and (g — 2) M constraints, see refs. [14]. We remark that most of the models which 
satisfy the WMAP bounds, the plus-shaped and squared points in the figures, satisfy 
the Higgs funnel condition tua — 2m x (see fig. 2, left panel) of ref. [10]): the main 
annihilation channel is XX ~> bb/rf through s-channel heavy Higgs exchange which is 
strongly constrained in our scenario and thus it is natural to have a reduction of photons 
emission. A reduction of the gamma ray flux in the funnel region it also found in ref. [3] 
in the study of the minimal supergravity (mSUGRA) plus right-handed neutrinos for 
see-saw generation of neutrino masses respect to the case of pure mSUGRA. 

4. Summary and conclusions 

In the framework of the MSSM with heavy SUSY-QCD particles and large tan/3 
we have studied lepton flavor violation in r — fi sector mediated by the heavy neutral 
Higgs $ = A, H at high energy colliders through the production and decay at LHC, 
pp — > $ + X, $ — > T[i and the \i — r fusion at a photon collider, 77 — > T[ibb. 

We have found that in models with 0.09 < n x h 2 < 0.13 and 1 < aff SSM x 10 9 < 4: 
at LHC the cross section for pp — > $ + X, <f> — > t/j, can reach O(10^ 1 — 10~ 2 ) fb in the 
range tua — 400 — 1000 GeV giving up to 10 events with 100 fb -1 ; the cross section of 
77 — > T/j,bb reaches O(10 -3 ) fb, thus too small even for the large value of the expected 
luminosity of 500 fb -1 . Prospects are somewhat more encouraging if we relax the lower 
limits, imposing only ft x h 2 < 0.13 and aff x 10 9 < 4: the cross section at LHC is 
about 2 fb for low tua masses and around 2 x 10~ 2 fb in 77 collisions. On the other 
hand, to observe such effects, in any case, the full luminosity of the machine is needed. 

We have also studied the spin-independent neutralino nucleus cross section: we have 
shown that in models that satisfy 0.09 < tt x h 2 < 0.13 and 1 < a^ SSM x 10 9 < 4, the 
cross section lies just below the sensitivity of XENON100 which should report results 
soon. The full XENON 1 ton is needed to cover all the parameter space. For models 
with a relic density inside the WMAP interval the particle physics factor for gamma rays 
flux from neuralino annihilation in the haloes is found to be smaller than the typical 
values found in similar studies without the constraints from lepton flavor violation and 
the updated B physics and (g — 2) M constraints. 

* * * 
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